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1.  INTRODUCTION 


The  FELS  (Free  Electron  Laser  Simulation)  program  was  conducted  at 
Battelle,  Columbus  Division,  using  a  modified  Nd  glass  laser  to  obtain 
experimental  data  on  the  interaction  of  a  FEL  with  composite  materials.  A 
simultaneous  effort  at  PSI  was  directed  toward  the  development  of  a  basis  for 
understanding  the  performance  of  these  materials  and  the  extrapolation  of  the 
data  to  the  full  scale  scenario.  The  conclusions  of  this  study  (Simons  et  al., 
1989)  suggested  continued  data  analysis  and  model  development  to  be  completed 
on  an  interim  study  (IFELS)  prior  to  continued  experimental  work  on  the 
advanced  FELS  (AFELS)  program.  This  report  documents  the  results  of  the  IFELS 
study  at  Physical  Sciences  Inc.  (PSI). 

Since  the  final  report  (Simons  et  al.,  1989)  on  the  original  FELS  study 
was  organized  phenomenologically,  this  report  was  organized  similarly  in  order 
that  we  may  report  the  progress  in  each  area.  Thus,  Section  2  is  organized  to 
illustrate  the  accomplishments  in  surface/subsurface  phenomena  and  Section  3 
discusses  the  progress  in  the  plume  absorption  and  scattering  issues. 
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2.  SURFACE/SUBSURFACE 


The  efforts  in  this  area  of  the  phenomenology  have  concentrated  on  the 
radiometric  data  on  black  materials  (Subsection  2.1),  the  modeling  of  the 
surface  temperature  and  mass  loss  in  the  translucent  materials  (Subsection  2.2) 
and  on  extending  the  equilibrium  thermochemistry  to  include  larger  molecules 
(Subsection  2.3).  Each  of  these  studies  has  important  implications  when 
extrapolating  known  results  to  large  time  and  large  laser  spots. 

2.1  Black  Materials 

Testing  of  materials  under  the  FELS  program  utilized  sufficiently  short 
test  times,  i.e.,  sufficiently  low  fluence  on  target,  that  the  residual  energy 
left  in  the  target  significantly  affected  the  measured  value  of  Q*.  Hence,  it 
was  necessary  to  understand  the  residual  energy  measurements  in  order  to  scale 
the  data  to  a  true  steady  state.  It  has  long  been  assumed  that  the  surface 
temperature  and  residual  energy  are  controlled  by  the  same  physics.  Analysis 
of  the  FELS  data  (Simons,  et  al.,  1989)  revealed  that  the  residual  energy  is 
ten  times  greater  than  that  which  would  be  consistent  with  the  radiometry.  It 
was  concluded  that  residual  energy  is  controlled  by  that  part  of  the  beam 
(<f>=0.1)  that  "leaks"  between  the  fibers  and  is  deposited  over  depth  L  (tens  of 
microns)  below  the  surface.  This  model  was  confirmed  with  model/data 
comparisons.  The  surface  temperature,  however,  is  controlled  by  that  part  of 
the  beam  ( 1—  4>=0 . 9 )  that  is  deposited  on  the  fiber  in  the  absorption  depth 
and  diffuses  into  the  fiber  with  fiber  diffusivity,  Vf.  Surface  temperature 
predictions  were  presented  but  the  radiometry  data  was  not  yet  reduced  for 


comparison.  Some  of  this  radiometry  data  was  reduced  in  the  IFELS  program  and 
is  utilized  below  to  confirm  the  physics  of  the  model. 

The  predicted  minimum  surface  temperature  just  before  the  last  pulse  of 
each  pulse  train  (30  to  40  pulses)  is  illustrated  in  Figure  1.  The  temperature 
increases  with  pulse  duration  Tp  because  there  is  more  time  for  the  energy  to 
diffuse  into  the  fiber.  The  temperature  increases  with  decreasing  pulse 
fluence  because  there  is  less  fiber  recession  to  overtake  the  diffusion  front. 
The  WLT-1  data  (2  ps  pulses,  36  pulse  trains)  support  the  fluence  dependence  of 
the  model  but  the  pulse  duration  dependence  is  yet  to  be  confirmed.  The 
predicted  temperatures  are  slightly  higher  than  actually  measured.  These 
predictions  may  be  reduced  by  reducing  the  assumed  fiber  diffusivity  from 
0.1  cm^/s  to  0.05  cm^/s.  Figure  2  illustrates  the  consistency  of  this  approach 
in  predicting  the  minimum  temperature  at  the  end  of  the  train.  Radiometry  data 
(Figures  3  and  4)  for  the  temperature  between  the  first  three  pulses  of  the 
first  two  trains  (Pulse  1,  2,  3  and  37,  38,  39)  also  support  the  predictions 
illustrated  in  Figure  2. 

As  short  pulse  radiometry  data  becomes  available,  we  will  compare  it  to 
our  predictions.  Vhile  it  is  not  our  intention  to  make  surface  temperature 
predictions  per  say,  it  is  our  goal  to  validate  the  pulse  fluence  and  pulse 
duration  dependence  of  our  model  against  all  observables  in  order  that  we  may 
more  confidently  predict  the  corresponding  dependence  of  Q*  on  these  laser 
parameters . 

2 . 2  Translucent  Materials 

During  the  original  FELS  program,  PSI  (Simons,  et  al.,  1989)  developed 
models  for  time  dependent  charring  and  explosive  pyrolysis  of  translucent 
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Figure  4. 

Front  Surface  Temperature  Between  Pulses  37  and  39 
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materials.  At  the  conclusion  of  the  FELS  program,  our  efforts  were  directed 
toward  combining  these  effects  into  a  single  computational  capability  for 
predicting  observables  such  as  mass  loss,  recession,  surface  temperature  and 
residual  energy.  However,  several  of  the  model  inputs  are  still  not  accurately 
known:  the  most  critical  are  the  fluence  dependent  absorptance  of  the 
material,  the  absorption  depth  of  the  fully  charred  material,  the 

pyrolysis  rate  (k)  of  the  resin,  the  temperature  dependent  strength  au(T)  of 
the  material,  and  the  enthalpy  of  the  pyrolysis  gas  when  exposed  to  the  surface 
temperature  Ts. 

Baseline  values  of  these  parameters  have  been  determined  and  sensitivity 
studies  carried  out  to  isolate  the  most  important  unknowns  in  order  that  we  may 
define  a  test  series  to  either  directly  or  indirectly  measure  them.  The 
baseline  values  of  these  parameters  are  discussed  below. 

The  enthalpy  of  the  pyrolysis  gas  when  exposed  to  the  surface  temperature 
Ts  and  corresponding  surface  pressure  have  been  calculated  assuming  equilibrium 
thermochemistry.  To  a  good  approximation,  the  pyrolysis  gas  behaves  as  though 
it  were  a  gas  with  a  Cp  of  approximately  5  up  to  a  temperature  of  4000  K  with  a 
constant  enthalpy  (-20  kJ/g)  above  4000  K. 

The  resin  pyrolysis  rates  were  obtained  by  extrapolating  measured  values 
(obtained  at  1000  K  and  heating  rates  of  1  K/s)  to  temperatures  of  5000  K  and 
heating  rates  of  10^  K/s.  These  extrapolations  are  illustrated  in  Figure  5 
and  are  supported  by  extrapolated  coal  (a  simmilar  C-0-H  system)  pyrolysis  data 
taken  at  2000  K  and  heating  rates  of  105  K/s.  A  factor  of  ten  uncertainty  in 
the  pyrolysis  rate  is  evident. 
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Figure  5. 

Extrapolated  Pyrolysis  Rates 
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PYROLYSIS  TIME 


The  material  strength  and  thermal  degradation  thereof  are  assumed  to  be  of 
the  form 

«r(T)  =  a  exp(-T/300)  exp(T/T) 

where  a0  (1000  atm)  is  the  yield  stress  at  room  temperature  and  T  =  1000  allows 
that  to  decay  to  60  atm  at  blowoff  temperatures  of  the  order  of  2000  K. 
Similarly,  T  =  1500  reduces  the  strength  to  15  atm  at  that  same  blowoff 
temperature. 

In  working  with  the  explosive  pyrolysis  model,  it  is  important  to  note 
that  the  gas  pressure  in  the  pores  increases  with  -fie.  When  the  gas  pressure 
exceeds  the  yield  stress  cru ,  mechanical  mass  loss  occurs.  Hence,  the  ratio  of 
4F  to  <ru  is  the  critical  parameter  for  establishing  mechanical  mass  removal. 

The  average  absorptance  is  yet  another  major  variable.  Measurements 
(Figure  6)  illustrate  that  the  absorptance  is  as  low  as  0.25  on  the  virgin 
material  and  as  high  as  0.7  on  the  fully  charred  material.  The  absorptance  is 
assumed  to  vary  linearly  with  the  fraction  (f)  of  the  material  charred, 
a  =  0.25  (1-f)  +  0.7  f 

Since  f  varies  with  depth,  we  must  determine  the  proper  average  value  of  a  to 
use.  Two  approaches  have  been  used:  1)  the  evaluation  of  f  at  the  material 
surface,  and  2)  the  weighting  a(f)  with  the  change  in  intensity  in-depth. 

1)  <x(o)  =  a  at  x  =  0 

2)  a=  /  adl/J  dl 

The  baseline  values  of  these  parameters  may  be  determined/confirmed 
through  comparison  of  the  model  to  FELS-7C  data  for  material  110. 

Figures  7  and  8  illustrate  the  effect  of  on  surface  temperature  and 
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Figure  6. 

Post-Irradiation  Optical  Properties 
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Figure  8. 


Determination  of  From  Mass  Loss  Data 


mass  loss  respectively.  Larger  values  of  X-y®  cause  more  blowoff  (lower  Q*) 
and  colder  surface  temperatures.  A  value  of  6  ym  is  consistent  with  data 
wherein  the  surface  temperature  is  the  more  sensitive  variable. 

The  material  strength  parameter  T  is  varied  to  illustrate  the  correspond¬ 
ing  changes  in  surface  temperature  and  mass  loss  in  Figures  9  and  10 
respectively.  In  these  cases,  blowoff  occurred  at  2000  K.  Hence  T  =  1000  and 
1500  K  correspond  to  ultimate  strengths  of  60  and  15  atm  respectively.  Surface 
temperature  is  again  the  more  sensitive  variable  and  suggests  a  value  of  T  of 
the  order  of  1000  K. 

The  role  of  pyrolysis  rate  is  illustrated  in  Figures  11  and  12.  Higher 
pyrolysis  rates  increase  the  blowoff  and  expose  colder  surfaces  while  lower 
pyrolysis  rates  reduce  blowoff  and  cause  the  surface  material  to  remain  in  tact 
and  be  heated  to  elevated  temperatures.  The  extrapolated  values  of  pyrolysis 
rates  are  consistent  with  both  the  available  data  and  the  current  values  of  T 
and  X^j®. 

All  of  the  calculations  presented  above  predict  a  sharp  increase  in  Q*  at 
fluences  above  10  J/cm^.  This  is  due  to  the  enthalpy  of  the  pyrolysis 
gas.  The  sensitivity  of  the  calculated  mass  loss  to  this  property  of  the 
target  material  is  illustrated  in  Figure  13.  Using  equilibrium  to  determine 
the  effective  value  of  Cp,  we  obtained  Cp  ~  5  whereas  for  a  non-reacting  gas, 

Cp  is  of  order  2.  The  mass  loss  data  support  the  non-reacting  limit  in  which 
the  pyrolysis  products  (large  aromatic  rings)  do  not  decompose  further  into  C, 
C2,  and  C3.  Future  test  series  which  measure  plume  gas  species  will  illuminate 
the  nature  of  the  pyrolysis  gas  and  the  corresponding  value  of  the  enthalpy. 
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Figure  13. 

Sensitivity  of  Mass  Loss  to  Pyrolysis  Gas  Enthalpy 
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The  sensitivity  of  the  surface  temperature  and  mass  loss  predictions  to 

the  material  absorptance  is  illustrated  in  Figures  14  and  15  respectively. 

Clearly  the  single  most  important  parameter  controlling  mass  loss  at  low 

fluence  (Fp  <  10  J/cm^)  is  the  relationship  between  charring  and  material 

absorptance.  A  complete  understanding  of  the  data  illustrated  in  Figure  6  will 

be  fundamental  to  the  understanding  of  this  phenomena. 

The  parameter  studies  described  above  have  indicated  that  the  time 

dependent  charring  and  absorptance  are  the  major  factors  controlling  mass  loss 

in  a  translucent  material  at  low  fluence.  The  absorptance  and  reflectance  of  a 

translucent  material  are  controlled  by  the  scattering  at  the  fiber/resin 

interface  and  a  two  flux  model  is  necessary  to  describe  this  phenomena. 

If  the  forward  flux  is  i  and  the  backseat tered  is  j,  then  the  i  component 

is  decreased  due  to  scattering  S  and  absorption  a  but  enhanced  due  to 

scattering  of  the  j  component.  Hence 

di  „ .  .  _ . 

Hx  =  "  Si  ~  “i  +  SJ 

and 

=  -  Sj  -  <xj  +  Si 

where  x  =  -y.  If  S  and  a  are  constant  then  i  and  j  are  of  the  form 
i  -  I0 

and 

j  =  R  IQ  e-fr 

where  the  effective  absorptivity  3  is 
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ce  Temperature  to  Absorptance 


and  the  reflection  coefficient  R  is 


R  = 


|3  -  a 

3  +  a 


Rewriting  the  relationship  between  S  and  R 


or 
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for  the  virgin  material  (Rv  =  0.75)  and 


2S 

c 


a 

c 


2.45 


for  the  charred  (Rc  =  0.3)  material.  The  corresponding  relationship  between  (3C 
and  a,,  is  then 

3C  =  1.86  a,, 

where  3c  is  identically  the  reciprocal  of  the  effective  absorption  depth  (6  ym) 
of  the  fully  charred  material  in  the  single  flux  model  used  in  the  parameter 
studies  above 

3C  =  (6  ym)-1 

and  Oc  now  reflects  the  actual  absorption  depth  of  the  fully  charred  material 


X  1  1-86 

v®  a  3 

c  c 


10  ym 


Similarly,  for  the  virgin  material, 
Pv  =  7  “v 
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and  the  virgin  absorptivity  is  approximately 
c<y  =  1/X-vjq  =  1  cm-1 

vhich  requires  an  effective  absorption  depth  of 
Ay  =  l/|3y  =  1.4  mm 

due  to  the  multiple  passes  in  the  scattering  environment. 

Assuming  that  scattering  occurs  on  the  interface  between  the  resin  and  the 
fiber,  we  write 

S  =  2cr/D 

where  we  have  assumed  that  there  are  two  surfaces  per  lattice  spacing  D 

(D  =  10  pm)  and  a  is  the  reflection  per  surface.  The  corresponding  properties 

of  the  virgin  and  charred  materials  are: 


Rv  =  0.75 

Rc  =  0.30 

ay  =  1  cm-3 

a c  =  103  cm-1 

Py.  =  7  cm-1 

3C  =  1860  cm'1 

Sv  =  24  cm-3 

Sc  =  1225  cm"1 

Oy  =  0.012 

ac  =  0.61 

Using  the  properties  of  the  virgin  and  charred  material,  it  is  reasonable 
to  assume  that  the  absorbtivity  a  and  reflectance  R  scale  in  direct  proportion 
to  the  fraction  (f)  of  material  that  has  charred 
a  =  a^f  +  ay(l-f) 

R  =  Rcf  +  Rv(l-f) 

where  f  is  controlled  by  the  growth  of  N  sites  of  radius  r 
f  =  (4it/3)  r3N 

and  each  site  heats  and  grows  due  to  the  absorption  of  laser  irradiation. 


We  previously  balanced  the  absorbed  intensity  1^  with  the  heating,  growth 


and  spherical  conduction  of  energy 

ml  h  *  It  (j  nr3  pcT)  +  kt  4nr2  ©  ' 

If  the  site  heats  (or  cools)  when  T  <  Tmax  and  grows  when  T  =  Tmax  we  may  write 
9T  rIA-4ktT 

at  "  ~ 5  2 

3  pcr 


for  r  fixed  and  T  <  Tmax  or 


»  rl.  -  rk  T 
3r  A  t  max 


at 


4r  pcT 


max 


for  T  -  Tmax  and  3r/3t  >  0.  Currently,  1^  is  expanded  to  include  both  the  i 
and  the  j  component 

IA  =  gdi  +  Ij) 

and  a  geometric  factor  g  to  account  for  the  fact  that  the  initial  site  (r  =  r0) 
is  completely  absorbing  and  the  charred  material  around  the  site  absorbs  with 
depth  i  « e  *  , 

jit^Ia  =  itr0^(Ii+Ij ) 

+  ntr^r^XIi+IjXl-e-r/Xvco) 

or 

g  =  1  -  (1  -  rQ2/r2)  G-r/Xv«> 

The  two  flux  deposition  model  has  been  programmed  into  the  PSI  ablation/ 
explosive  pyrolysis  code.  Initial  de-bugging  is  complete  and  the  model  is 
undergoing  initial  model/data  comparison  with  the  FELS-7B,  7C  and  WLT-1  mass 
loss  and  radiometry  data. 
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2.3  Extension  of  Equilibrium  Calculations 


2.3.1  Introduction 

Equilibrium  thermochemistry  has  been  used  to  predict  gas  species  during 
vaporization  and  the  subsequent  values  of  Q*.  However,  the  calculations  did 
not  include  (for  lack,  of  thermodynamic  data)  large  carbon  rings  as  an  allowable 
species.  The  presence  of  large  amounts  of  aromatic  rings  would  not  only  reduce 
Q*,  but  would  impact  the  work  on  beam  absorption  by  gas  species  as  well. 

Below,  we  extend  our  previous  analysis  to  include  the  effect  of  Cg,  C^q,  cgH, 
CqH,  and  CiqH  on  the  equilibrium  properties  of  graphite  epoxy  and  epoxy  under 
pyrolysis  and  vaporization  conditions.  First  we  shall  examine  the  thermo¬ 
dynamic  properties  of  these  large  molecules  and  describe  our  method  of  generat¬ 
ing  approximate  partition  functions  for  them.  Comparison  is  made  with  existent 
Sandia  data.  Secondly,  we  will  examine  under  what  conditions  of  temperature, 
pressure,  and  elemental  composition  these  molecules  contribute  significantly  to 
the  equilibrium.  The  constrained  equilibrium  case  is  also  examined.  The  rela¬ 
tive  importance  of  CgH  to  Cg,  CgH  to  Cg  and  C^qH  to  C^q  is  illustrated.  It 
should  be  pointed  out  that  the  spirit  of  this  work  was  exploratory  in  nature. 

As  a  result  some  of  the  approximations  made  can  be  readily  and  obviously 
improved.  These  potential  improvements  will  be  pointed  out  along  the  way. 

2.3.2  Approximation  to  Free  Energy  and  Formation  Enthalpy 

To  generate  curve  fits  for  the  partition  functions  for  the  molecules  Cg, 
C10,  CgH,  CjqH,  and  CgH,  information  on  free  energy  in  the  form 
[G°(T)  -  H°(Tr)]/T  and  formation  enthalpies  at  the  reference  temperature 
Tr  =  298.15  K  are  required.  Lacking  thermodynamic  data  we  make  the 
approximations  to  both  the  free  energy  and  formation  enthalpy: 
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1. 


Linearly  extrapolate  from  known  data  (Gustafson,  1986;  Lee  and 
Sanborn,  1973;  Leider  et  al.,  1973;  Chase  et  al.,  1985)  for  to 

larger  C-clusters. 

2.  Treat  CgH  as  Cg  and  CiqH  as  C^o>  i.e.,  neglect  the  modifications 
introduced  by  the  C-H  bond  entirely  except  in  the  composition  and 
molecular  weight;  the  energetics,  however,  remain  unchanged. 

The  extrapolated  free  energy  for  seven  different  temperatures  appears  in 
Figure  16,  while  the  extrapolated  formation  enthalpy  at  298  K  appears  in 
Figure  17.  The  upper  limit  of  6000  K  conforms  to  well  established  information. 
As  an  example,  the  extrapolated  values  for  C^q  appear  in  Table  1  along  with 
thermochemical  information  and  the  curve  fit  for  ln(Q)  as  a  function  of  temper¬ 
ature  generated  by  a  non-linear  least  squares  fitting  routine.  The  chemical 
equilibrium  code  used  requires  enthalpies  relative  to  absolute  zero  degrees  not 
room  temperature  T  =  298.  This  enthalpy  increment  which  is  10.51  kJ  mol-^-  and 
amounts  to  about  0.5  percent  of  the  formation  enthalpy  of  C^q  at  298  K  has  been 
neglected.  As  can  be  seen  the  H  and  Cp  fits  have  unphysical  negative  values. 
This  is,  however,  not  surprising  since  we  have  used  seven  values  to  generate  a 
seven  parameter  fit.  An  obvious  improvement  can  be  made  here. 

2.3.3  Comparison  with  Existing  Information 

We  compare  the  inputs  to  the  partition  function  fitting  routine,  namely 
-[G°(T)  -  H°(298)]/T  and  AHf°(298)  to  Sandia  thermochemical  fits  (Burcat,  1984; 
Miller  and  Bowman,  1988).  The  results  of  our  comparison  appear  in  Tables  2 
and  3.  Comparing  CgH  we  find  a  40  percent  discrepancy  in  the  formation 
enthalpy  but  the  free  energies  are  surprisingly  close  (Table  2).  For  CgH  we 
find  an  18  percent  discrepancy  in  the  formation  enthalpy  and  a  7.5  percent 


27 


1  -  600  °K 
2-1000 

3  -  2000 

4  -  3000 

5  -  4000 

6  -  5000 

7  -  6000 


4  5  6  7  8  9 

NUMBER  OF  ATOMS  N 

B-2747 


Figure  16. 

Free  Energy  for  Ci~C;]o 


Table  1.  Extrapolated  Thermo chemical  Properties  of  C3g 


Molecule:  C3g,  mw  =  120.11,  Ideal  Gas 

AHf (298. 15)  =  1600.0  KJ  mol"1 

T 

[G°(T)-H°(298)]/T 

Cp  Fit 

[H°(T)-H°(298)]Fit 

AGf°(T) 

(K) 

JK"1  mol'1 

In  Q  Fit 

JK-1  mol"1 

KJ  mol"1 

KJ"1  mol-1 

600.0 

450.0 

-2.22654E+02 

-1.96492E+02 

-4. 18624E+02 

1377.76 

1000.0 

510.0 

-8.76575E+01 

2.48971E+02 

-3.69541E+02 

1216.62 

2000.0 

620.0 

2. 10979E+01 

1.65496E+02 

-1.65293E+02 

820.16 

3000.0 

690.0 

6. 11844E+01 

3.06264E+02 

5.48397E+01 

453.31 

4000.0 

750.0 

8.41495E+01 

4 . 15464E+02 

4.37877E+02 

76.52 

5000.0 

805.0 

1.00163E+02 

7.44302E+00 

7 . 13547E+02 

-325.15 

6000.0 

840.0 

1 . 10605E+02 

-1.42778E+03 

1 . 10408E+02 

-658.70 

Partition  Function  Fit  In 

LX  =  -1.83308E+03 

L2  =  2.71180E+02 

L3  =  -2.03945E+04 

L4  =  -3.02675E+07 

L5  =  -1.32798E-01 

L6  =  1.57129E-05 

L 7  =  -8.65216E-10 

Q  =  Li  +  L2  InT  -  L3/T  +  L4/T2  +  L5T  +  L5T2  +  L7T3 

I  1 


(i  ii  ii  it  (a  ii  is  ii  ii  ii  fi  ia  ii  ii  i  i  f«  ii  ia 


Table  2. 


Comparison  of  PSI  and  Sandia  Thermochemical  Data 


Molecule:  CgH 


AHf (298)  kJ  mol-1 

PSI  Sandia 

1255.2  891.93 


PSI 

Sandia 

T 

[G°(T)-H°(298)]/T 

T 

[G°(T)-H°(298)]/T 

(K) 

JK"1  mol-1 

(K) 

JK"1  mol'1 

200 

3. 16600E+02 

298 

310.102 

298 

3.08780E+02 

300 

310.104 

400 

3. 12420E+02 

400 

313.899 

600 

3 . 27650E+02 

500 

321.473 

800 

3.44510E+02 

600 

330.391 

1000 

3.60540E+02 

700 

339.707 

1200 

3. 75350E+02 

800 

349.003 

1400 

3.88990E+02 

900 

358.085 

1600 

4.01580E+02 

1000 

366.866 

1800 

4. 13210E+02 

1500 

405.910 

2000 

4.24010E+02 

2000 

438.110 

2200 

4 . 34090E+02 

2500 

465.306 

2400 

4. 43500E+02 

3000 

488.772 

2600 

4.52330E+02 

3500 

509.364 

2800 

4 . 60620E+02 

4000 

527.679 

3000 

4.68480E+02 

4500 

544.157 

3200 

4.75930E+02 

3400 

4.82960E+02 

3600 

4.89700E+02 

3800 

4 . 96100E+02 

4000 

5.02210E+02 

4200 

5.08060E+02 

4400 

5. 13670E+02 

4600 

5. 19070E+02 
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Table  3. 


Comparison  of  PSI  and  Sandia  Thermochemical  Data 


Molecule:  CgH 


AHf (298)  kJ  mol"1 


PSI  Sandia 

1425.02  1208.72 


PSI 

Sandia 

T 

[G°(T)-H°( 298) ]/T 

T 

[G°(T)-H°(298)]/T 

<K) 

JK"1  mol-1 

(K) 

JK"1  mol"1 

600 

3.80000E+02 

298 

328.053 

1000 

4. 30000E+02 

300 

328.055 

2000 

5.20000E+02 

400 

332.744 

3000 

5.80000E+02 

500 

342.186 

4000 

6.30000E+02 

600 

353.393 

5000 

6. 70000E+02 

700 

365.173 

6000 

7. 10000E+02 

800 

376.984 

900 

388.565 

1000 

399.795 

1500 

450.001 
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discrepancy  in  the  free  energy  (Table  3).  We  expect  our  extrapolated  results 
to  become  better  as  we  go  to  higher  carbon  clusters.  A  better  comparison  to 
the  Sandia  fits  would  be  plots  of  InQ,  Cp,  H-H(298),  and  AGf  as  a  function  of 
temperature  which  are  readily  generated. 

2,3.4  When  are  the  Large  Molecular  Clusters  Important? 

We  examine  under  what  conditions  of  T,  P,  and  composition  large  carbon 
clusters  are  important.  For  definiteness  we  focus  on  but  other  clusters 
have  the  same  qualitative  features.  We  then  examine  how  CiqH  fares  against 
C^Q'  We  shall  find  that,  most  of  the  time  «  1,  and  hence  is 

usually  not  an  important  species  under  equilibrium  conditions.  The  low  tem¬ 
perature  equilibrium  is  dominated  by  the  largest  cluster,  while  at  higher 
temperatures  there  is  a  competition  between  the  various  carbon  molecular 
species.  We  shall  see  we  can  stabilize  to  some  extent  the  heavier  molecules  by 
constraining  the  equilibrium  and/or  varying  the  composition. 

The  temperature  dependence  of  the  mole  fraction  of  C^q  at  1  atm  is  shown 

in  Figure  18  (solid  curve)  for  graphite  epoxy  under  pyrolysis  conditions.  [It 

is  of  interest  to  make  similar  plots  for  increasing  cluster  size  to  determine 

the  temperature  at  which  the  most  stable  cluster  disappears.]  Above  4000  K 

there  is  virtually  no  C^o  at  equilibrium,  while  below  this  temperature,  up  to 

3100  K,  it  is  the  dominant  pure  carbon  molecular  species.  Above  3100  K  C3 

becomes  more  dominant  and  eventually  other  molecular  carbon  species  compete 

(see  Figure  19).  In  Figure  20  the  enthalpy  is  displayed:  H  (T=500)  = 

L10 

14.35  kJ/g,  Hr  (T=3000)  =  17.96  kJ/g,  which  is  substantially  higher  than  what 
C10 

one  would  get  if  solid  carbon  were  condensed  (the  stable  phase  under  these 


MOLE  FRACTION 


TEMPERATURE  (K) 


Figure  18. 

Mole  Fraction  of  C^q  at  P  =  1  atm  in  Graphite  Epoxy  under  Equilibrium 
(Solid  Curve)  and  Constrained  Equilibrium  (Dashed  Curve) 
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B-2751 


Figure  19. 

Graphite  Epoxy  Equilibrium  Composition  (1  atm) 
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conditions)  as  is  expected.  Above  4000  K  the  enthalpy  agrees  with  that  of  a 
previous  report  (Resendes,  1987). 

Let  us  now  consider  a  fixed  amount  of  and  study  the  pressure  depen¬ 
dence  of  the  vapor  as  we  vary  the  temperature.  We  examine  the  equilibrium  (any 
other  equilibrium  involving  C^q  fine): 


iOC  *  C10 


(1) 


along  with  the  consequent  law  of  mass  action  which  we  write  in  the  form: 


TO  P 
TO"  =  "9 


u°_  (T)  -  10u°  (T) 

TO  C 

RT 


(2) 


In  Eq.  (2)  P  denotes  the  total  pressure  in  atms,  subscripted  u°  the  molar  free 

energy  and  subscripted  n  the  number  of  moles.  Provided  the  total  number  of 

moles,  n-p,  does  not  change  "much"  with  P  and  T  any  other  equilibrium  at  a 

different  P  and  T  such  that  n_  /n  =  constant,  necessarily  has  the  same  mole 

C10 

numbers.  For  the  sake  of  calculation  let  us  take  T  =  3000  and  P  =  1  atm  where 

the  mole  fraction  of  C^n  is  approximately  20  percent,  or  n  /nr  =  632.3629. 

TO 

Requiring  approximately  the  same  n_  /n  ratio  at  4000  and  5000  K,  we  find  the 

TO  C 

corresponding  pressures  using  the  previous  equation.  Note  nj  =  0.218442  (from 
code  output  at  3000  K)  and  R  =  8.3144  x  10“^  kJ/(mole  K).  Our  results  are 
plotted  in  Figure  21.  Comparison  of  these  number  with  the  corresponding 
vaporization  conditions  of  Ref.  3  shows  that  for  3000  <  T  <  5000  there  is  a 
negligible  mole  fraction  of  C^q  under  vaporization  conditions.  It  also 
indicates  that  if  we  are  to  see  a  significant  amount  of  C^q  ( s20  percent  mole 
fraction)  at  these  temperatures  the  pressures  are  extremely  high. 
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TEMPERATURE  (K) 

B-2753 


Figure  21. 

Pressure  Dependence  on  Temperature  for  20  Percent  Mole  Fraction  of  C^q 
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To  see  the  effect  of  constrained  equilibrium  on  the  presence  of  large 
clusters  we  refer  to  Eq.  (2)  and  note  that  the  total  number  of  moles  n^  appears 


in  the  denominator  raised  to  the  nineth  power.  This  indicates  that  if  we 

constrain  the  equilibrium  by  not  allowing  the  smaller  species  to  appear  this 

will  reduce  n^  thereby  increasing  n  /nr.  To  illustrate  the  effect  we  will 

10 

not  allow  C,  CH,  CH2,  CO,  OH,  HCO,  C2H.  The  implication  is  that  at  the  begin¬ 
ning  of  the  laser  pulse  when  the  smaller  species  are  practically  non-existent 
due  to  kinetic  limitations,  there  should  be  more  C^q  present  than  an  equilib¬ 
rium  calculation  indicates.  Figure  18  (dashed  curve)  gives  the  results  of  such 
a  calculation.  It  is  interesting  to  note  the  lack  of  C^q  (in  favor  of  CH4)  at 
1000  K.  It  is  not  clear  whether  this  represents  a  real  physical  effect; 
thermodynamically  CH4  is  favored  over  C^q.  Note  for  example  that  at  3000  K  the 
constrained  run  has  an  almost  30  percent  increment  of  C^q  relative  to  the 
non-constrained  run. 

We  want  to  compare  the  relative  abundances  of  C10H  and  ^10*  Our  starting 
point  is  the  law  of  mass  balance  for  the  equilibrium  reaction: 


C10  +  H  -  C10H 


(3) 


which  we  write  in  the  form: 


C10H 


'10 


=  P  e 


UH 

RT 


lnTj 


(4) 


[We  are  considering  here  a  non-constrained  equilibrium. ]  Since  our  normaliza¬ 
tion  is  such  that  n-p  <  1,  ng/n^  <  1  always.  ug°  >  0  at  low  T  when  n^  is  small 
and  U[j  <  0  at  high  T  when  ng  is  large  with  the  result  that  (at  constant 
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pressure)  nr  g/nr  «  1 •  The  remaining  molecules  C5H  and  Cg,  CgH  and  Cg 
L10  C10 

follow  the  same  trend. 

Lastly,  we  consider  the  effect  of  modifying  the  C/H/O  ratios  on  the 
previous  results.  We  do  this  by  working  with  an  epoxy  rather  than  a  graphite 
epoxy.  From  Eq.  (4),  the  n^^/nQ  ratio  will  be  favored  in  epoxy  versus 
graphite  epoxy  if  nj  for  the  same  P,  T  and  equilibrium  species  is  lowered  in  the 
epoxy  relative  to  the  graphite  epoxy.  A  code  calculation  shows  that  n^  in  epoxy 
is,  in  fact,  higher  than  in  graphite  epoxy  so  that  ncjQ<epoxy)/nQjQ(g.e. )  <  1. 

However,  nu  is  higher  in  epoxy  so  that  n  u/n_,  (epoxy)  >  n„  u/n„  (g.e.). 

H  C1qH  C1Q  C1qH  C1Q 

2.3.5  Conclusion 

In  general,  what  we  find  is  that  at  a  given  pressure  the  higher  molecular  C 
species  are  favored  at  low  temperatures.  If  the  pressure  is  increased  they  can 
be  condensed  out.  Changing  the  composition  from  carbon  fiber  poor  to  fiber  rich 
favors  high  molecular  carbon  formation.  These  conclusions  present  no  surprises. 
Interesting  effects  arise  from  constraining  the  equilibrium  which  is  a  limited 
form  of  kinetics  modeling.  It  was  found  that  the  presence  of  high  molecular 
carbons  could  be  enhanced  by  disallowing  the  formation  of  smaller  chemical  spe¬ 
cies.  However,  one  can  find  regions  (in  P,  T  space)  where  this  is  not  the  case 
(recall  Figure  18).  Under  vaporization  conditions  the  amount  of  Ciq  present  is 
negligible.  Also  the  ratio  CnH/Cn  for  large  n  is  always  much  less  than  one.  All 
this  leads  to  the  conclusion  that  knowledge  of  the  chemical  distribution  as  well 
as  P,  T,  time  histories  are  necessary  for  a  proper  calculation  of  the  chemical 
contribution  to  Q*.  It  is  also  clear  that  extrapolation  of  the  type  that  was 
done  here  for  higher  molecular  clusters  will  yield  the  same  qualitative  results. 
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3.  PLUME  DEBRIS/BEAM  INTERACTION 


The  interaction  of  the  beam  with  gas  and  particulate  phase  debris  is 
potentially  important  in  scaling  Q*  with  spot  size.  The  particulates  are 
studied  in  detail  in  Subsection  3.1  and  are  potentially  capable  of  scattering 
the  beam  and  increasing  Q*  by  an  order  of  magnitude.  Similarly,  gas  phase 
absorption  studies  are  discussed  in  Subsection  3.2  in  order  to  assess  Q* 
increases  and  plasma  ignition  with  spot  size. 

3. 1  Scattering  by  Particulates 

The  previous  FELS  program  identified  that  a  few  percent  of  the  mass 
removed  from  a  target  was  present  in  submicron  particles  which  scattered  both 
the  1.06  pm  beam  and  the  6328A  HeNe  probe  beam.  The  HeNe  probe  beam 
experienced  attenuation  at  intensities  as  low  as  100  kW/cm^  whereas  the  Q* 
measurements  suggested  that  scattering  is  not  affecting  the  transmission  of  the 
1.06  pm  beam  to  the  target  at  intensities  below  1  MW/cm^.  An  analysis  of 
single  and  multiple  scattering  was  initiated  to  resolve  this  issue.  The 
updated  results  and  conclusions  are  reported  below. 

The  most  important  single  scattering  parameter  is  the  opacity  T  of  the 
particulates  in  the  plume  flowfield.  The  opacity  is  calculated  from  the  total 
column  density  M^y  of  mass  above  the  target 
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where  Is  is  the  average  intensity  on  the  spot  diameter  Ds,  V  is  the  velocity  of 
the  particles  (10^  cm/s),  and  Qs  is  the  value  of  Q*  in  the  absence  of  plume 
absorption  or  scattering  (30  kj/g).  If  mass  fraction  f  of  Mcy  is  broken  into  N 
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resin  fragments  of  density  p  (1.25  g/cc)  and  radius  a,  the  optical  opacity  (r) 
becomes 


T  =  No 

s 


3f  I  D 
s  s 

4  p  V  Q* 

s 


f 


where  as  is  the  optical  cross  section  of  a  particle  of  radius  a.  The  optical 
cross  sections  are  obtained  from  Mie  code  calculations  and  are  illustrated  in 
Figure  22.  The  cross  sections  are  functions  of  the  particle  size  and  laser 
wavelength.  For  submicron  particles,  the  cross  sections  (hence  the  opacities) 
for  the  HeNe  are  three  times  that  of  the  1.06  pm  beam. 

THeNe  =  3  T1.06 

This  important  distinction  between  the  HeNe  and  1.06  beams  was  overlooked  in 
the  FELS  program  and  is  rectified  below. 

Rewriting  THejqe  as 

THeNe  =  ^max  f  ^s  ^s  » 
the  maximum  attenuation  occurs  when 
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max 


3( a  / na^) 

4pVQ* 
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3  x  10  4  cm/V 


corresponding  to 


a 

s 


14  pm 


for  submicron  (a  =  0.3  pm)  particles. 
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Figure  22. 

Cross  Section  per  Unit  Volume 


The  HeNe  detector  has  a  narrow  field  of  view  and  any  photons  scattered  by 
particles  will  never  be  scattered  back  on  to  the  detector.  The  HeNe 
transmission  TjjeNe  is  therefore  statistically  related  to  the  opacity  by 
THeNe  =  exP(-THeNe)  • 

A  typical  HeNe  transmission  plot  is  illustrated  in  Figure  23  for  a  CW  pulse  at 
124  kW/cm^  on  material  130  using  a  0.8  cm  diameter  spot.  A  HeNe  transmission 
of  90  percent  is  consistent  with  the  above  analysis  if  /  =  0.03,  i.e.,  if 
3  percent  of  the  mass  is  in  submicron  particulates.  While  submicron  particles 
have  been  observed,  no  quantitative  measurement  supports  this  requirement.  It 
is  worth  noting  that  an  8  to  10  percent  mass  loading  is  required  to  explain  the 
HeNe  data  for  material  110. 

While  these  same  particles  are  scattering  the  1.06  pm  beam,  the  target 
will  have  a  larger  "collection  angle"  than  did  the  HeNe  detector.  A  two  flux 
(forward  and  reverse)  model  for  multiple  scattering  (Simons  et  al.,  1989) 
yields  the  beam  transmission  Tg 


_1 _ 

T1.06 


(5) 


where  a  is  the  backscatter  coefficient  (fraction  of  intensity  I  backscattered 
per  scattering  event).  Since  the  onset  of  the  HeNe  attenuation  occurs  at  an 

JL 

order  of  magnitude  lower  intensity  than  does  the  onset  of  the  beam  (Q  ) 
attenuation,  it  must  follow  that 


CT  T1.06  *  TO  ('CHeNe) 


for  all  intensities  Is.  When  the  distinction  between  Ti.06  an<*  THeNe  was  not 
considered,  we  were  led  to  the  erroneous  conclusion  that  o  =  0.1.  However, 
with  the  knowledge  that  for  submicron  particles, 
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Figure  23. 

HeNe  Tranmission  for  Material  130  at  124  kW/cm2 
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THeNe  =  3  T1.06 


it  follows  that  the  required  value  of  a  is  cr  =  0.3,  and  is  confirmed  below 
using  a  Monte  Carlo  scattering  code. 

The  Monte  Carlo  scattering  code  determines  the  fractional  penetration  of 
laser  light  through  a  cloud  of  scattering  particles.  The  code  is  one¬ 
dimensional  inasmuch  as  no  variation  in  the  particle  density  of  the  cloud  is 
assumed  in  the  directions  parallel  to  the  surface.  Given  this  fact,  the 
variation  normal  to  the  surface  is  also  irrelevant;  only  the  total  density  of 
scattering  particles  per  unit  area  is  significant  (e.g.,  if  the  cloud  is  twice 
as  thick  but  half  as  dense,  the  rays  will  travel  twice  as  far  between 
collisions,  but  the  net  penetration  will  be  the  same).  The  rays  are  traced 
using  three-dimensional  ray  tracing  techniques,  converting  to  multi-dimensions 
would  not  be  a  problem. 

The  code  works  by  tracing  probabilistic  paths  of  many  light  rays  and 
noting  the  fraction  that  penetrate  the  particle  cloud.  The  incident  laser 
radiation  is  taken  to  be  normal  to  the  surface.  A  ray  is  assumed  to  move  a 
randomly  picked  distance  chosen  such  that  a  large  number  of  numerical  rays  will 
statistically  have  the  same  falloff  with  distance  as  would  the  laser  light. 
Since  the  unscattered  radiation  will  fall  off  exponentially  as  exp(-x/mfp),  the 
travel  distance  is  chosen  to  be  distance  =  -mfp*ln(1.0  -  ran),  where  ran  is  a 
random  number  between  0  and  1  and  mfp  is  the  mean  free  path. 

Upon  collision,  the  ray  is  scattered  or  absorbed.  The  fraction  of  rays 
absorbed  will  be  Ea/(ES  +  Ea),  and  is  played  for  randomly.  Here  Es  is  the 
macroscopic  scattering  cross  section  and  Ea  is  the  macroscopic  absorption  cross 
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section.  The  scattered  rays  are  redirected  randomly  in  a  manner  consistent 
with  the  differential  scattering  cross  sections  of  the  scattering  particles. 

The  cross  sections  are  determined  by  calling  a  Mie  scattering  subroutine. 
For  a  given  size  particle  with  given  real  and  imaginary  refractive  indices,  the 
routine  determines  the  absorption  and  scattering  cross  sections  and  the 
probabilities  of  different  scattering  angles.  For  a  situation  in  which  the 
cloud  consists  of  particles  of  many  sizes,  the  routine  is  called  a  number  of 
times  to  compute  the  cross  sections  for  particles  in  a  number  of  size  bins,  and 
the  overall  absorption,  scattering,  and  differential  scattering  cross  sections 
per  unit  volume  are  computed. 

A  ray  is  traced  until  it  is  absorbed  or  leaves  the  cloud,  either  by 
penetrating  to  the  surface  or  by  scattering  away.  The  ray  is  then  scored  as 
one  of  four  categories: 

1)  Absorbed  by  the  cloud, 

2)  Scattered  away  by  the  cloud, 

3)  Penetrating  the  cloud  without  scattering,  or 

4)  Scattered  but  eventually  penetrating  the  cloud. 

The  fraction  of  light  in  each  of  these  categories  is  computed  by  dividing  the 
rays  in  each  bin  by  the  total  number  of  rays  traced.  A  significant  number  of 
rays  must  be  traced  to  give  good  statistics.  The  uncertainties  scale  as  the 
square  root  of  the  number  of  rays.  For  example,  if  1000  rays  are  traced  and 
10  percent  penetrate  (100  rays),  then  the  standard  deviation  is  10  rays  and  the 
fraction  penetration  is  between  9  and  11  percent  to  one  standard  deviation.  No 
variance  reduction  procedure  is  used  since  for  a  simple  problem  like  this  the 
computation  times  are  modest. 
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The  beam  transmission  TB  versus  the  plume  opacity  fi.06  is  illustrated  in 
Figure  24.  The  effect  of  particle  size  dominants  the  opacity  T^.06’  ^ut 
otherwise  does  not  strongly  influence  the  tranmission.  The  Monte  Carlo  code 
clearly  supports  the  two  flux  model  (Simons  et  al.  ,  1989)  for  a  =  0.3. 

Using 

Is  =  I  Tg  , 

T1.06  may  also  be  rewritten  as 

*1.06  =  K0  /  Ds  I  TB  (6) 

where  K0  for  submicron  particles  is 
K0  =  10~4  W/cm 

such  that  remains  1/3  of  xBeNe. 

Solving  Eqs.  (5)  and  (6)  for  the  beam  transmission  in  the  large  attenua¬ 
tion  limit,  we  obtain 

Tb  =  FwTs1  ' 

Since  Q*  is  QS*/TB,  we  illustrate  the  role  of  scattering  through  the  enchance¬ 
ment  of  Q*/Q*.  A  typical  composite  must  produce  3  percent  (/  =  0.03)  of  the 
mass  loss  in  submicron  fragments  in  order  to  explain  the  HeNe  transmission 
data.  Figure  25  illustrates  that  for  /  =  0.03  in  the  laboratory  (Ds  =  1  cm), 
intensities  of  a  few  megawatts  are  required  to  observe  an  increase  in  Q*  due  to 
scattering.  If  the  composite  were  constructed  with  30  percent  of  the  mass  in 
submicron  scatterers,  then  laboratory  measurements  (Ds  =  1  cm,  /  =  0.30)  would 
observe  a  significant  Q*  enhancement  at  intensities  of  one  megawatt.  Similar 
materials  in  the  robust  SDI0  scenario  (Ds  =  10  to  100  cm)  would  experience  an 
order  of  magnitude  enhancement  in  Q*.  Hence,  the  issue  of  scattering  has  an 
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Tb  -  beam  transmission 


Figure  24. 

Scattering  Cloud  Transmissivity 
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I  -  AVERAGE  INTENSITY  (W/cm2) 


Figure  25. 

Q*  Enhancement  via  Scattering 


important  impact  on  the  use  of  lasers  in  the  robust  SDIO  scenario  and  labora¬ 
tory  data  at  a  few  MW/cm^  is  required  to  further  support  these  conclusions. 

The  above  analysis  was  carried  out  for  a  CW  plume.  Plumes  generated  in 
the  laboratory  may  be  distinctly  different  in  the  RP  and  CW  modes.  For  suffi¬ 
ciently  short  pulses,  debris  leaves  the  surface  in  a  "shell"  with  velocity  V 
(V  ~  10^  cm/s).  The  di  stance  the  shell  travels  between  pulses  at  rep  rate  f  is 
V/f  and  the  number  of  shells  in  spot  size  Ds  is  Dsf/V.  Hence,  when  Dsf/V  >>  1, 
the  rep  pulse  plume  will  look  like  a  CV  plume. 

Typical  laboratory  experiments  are  conducted  on  a  1  cm  spot  at  10  kHz 
(i.e.,  Dsf/V  =  1)  and  will  clearly  demonstrate  a  rep  pulse  behavior.  This  is 
demonstrated  with  a  HeNe  probe  beam  transmission  plot  illustrated  in  Figure  26. 
The  transmission  is  unity  in  the  absence  of  any  debris  plume.  At  t  =  0,  a 
single  pulse  reduces  the  transmission  to  zero.  The  transmission  recovers  in 
time  as  the  flow  expands  and  achieves  a  value  of  0.8  when  a  second  pulse  again 
reduces  it  to  zero.  The  process  continues  until  a  repetitive,  periodic 
flowfield  is  established  where  transmission  varies  between  some  minimum  value 


just  after  the  pulse  to  some  maximum  value  that  will  be  seen  by  the  pulse. 

The  time  evolution  of  the  opacity  generated  in  the  rep  pulsed  mode  at 
frequency  f  may  be  expressed  (Simons,  et  al.,  1989)  as 


where  M  +  1  is  the  total  number  of  laser  pulses,  each  of  which  generates  an 
identical  shell  of  debris  with  initial  opacity  T0.  This  expression 
satisfactorily  explains  the  typical  transmission  versus  time  relationship 
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observed  experimentally  and  is  used  below  to  assess  the  similarity  between  RP 
and  CW  plumes  on  large  laser  spots,  i.e.,  when  Dsf/V  >>  1. 

We  wish  to  compare  the  various  RP  modes  at  fixed  average  intensity. 

Hence,  when  f  is  increased  by  a  factor  of  three  the  corresponding  pulse  fluence 
and  corresponding  initial  opacity  xQ  must  be  reduced  by  a  factor  of  three. 
Holding  Ds  and  V  fixed  at  1  cm  and  104  cm/s  respectively,  we  compare  predicted 
HeNe  probe  beam  transmissions  in  Figures  27  to  29  for  rep  rates  of  10,  30,  and 
90  kHz  respectively.  The  corresponding  values  of  xQ  were  chosen  to  be  2, 
0.66667  and  0.22222.  These  calculations  clearly  demonstrate  the  RP  behavior 
when  Dsf/V  =  1  and  a  CW  like  behavior  for  Dsf/V  =  9.  Thus,  the  predictions  and 
conclusions  associated  with  Figure  25  above  apply  to  a  RP  mode  as  well  as  a  CW 
mode  provided  the  spot  size  and  rep  rate  satisfy  Dsf/V  >  10. 

3 . 2  Gas  Heating  from  Laser  Absorption 
3.2.1  Introduction 

In  the  previous  report  (Simons,  et  al.,  1989,  Subsection  3.3)  we  pointed 
out  the  importance  of  gas  heating  when  the  ablation  plume  absorbs  the  laser 
energy,  as  well  as  showing  that  the  rapid  ablation  of  target  material  at  high 
irradiances  will  produce  a  large  density  of  ablation  products  in  the  plume  in 
front  of  the  target.  In  this  high  density  plume,  gas  absorption  will  increase 
Q*,  the  energy  to  ablate  one  gram  of  material.  The  most  important  consequence 
of  the  plume  heating  is  that  the  increased  temperature  gives  a  higher  electron 
concentration,  and  the  higher  electron  concentration  gives  increased  absorption 
via  the  inverse  bremsstrahlung  processes,  which  leads  to  a  thermal  runaway  that 
produces  a  plasma.  As  plasma  production  causes  a  large  increase  in  Q*,  we  see 
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Figure  28. 

Probe  Beam  Transmission  for  Dsf/V  =  3 


B-2824 
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Figure  29. 

Probe  Beam  Transmission  for  Dsf/V  =  9 
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that  it  is  important  to  understand  when  plume  absorption  leads  to  plasma 
production. 

A  major  issue  is  whether  the  ablated  species  absorb  the  laser  light.  That 
is,  do  they  have  states  with  excitation  energies  corresponding  to  the  energy  of 
the  laser  photons.  Multiple-photon  effects,  where  multiple  photons  are 
absorbed  at  the  same  time,  will  be  important  when  the  I  values  are  large 
enough,  but  we  ignore  these  and  only  consider  single-photon  processes.  Atomic 
species  will  only  have  a  small  number  of  transitions  from  low-lying  states  at 
wavelengths  near  1.06  pm,  and  thus  are  less  likely  to  absorb  than  are  molecular 
species,  by  which  we  mean  both  stable  molecules  and  radicals.  In  the  0.4  to 
1.06  pm  region,  absorption  can  excite  electronic,  vibrational,  and/or 
rotational  states.  Most  electronic  states  have  higer  energies  that  require 
radiation  with  shorter  wavelengths  than  1.06  ym,  and  the  large  absorption  cross 
sections  for  vibrational  transitions  are  at  longer  wavelengths.  The  absorption 
lines  have  finite  widths,  and  high  pressure  will  broaden  the  lines.  Absorption 
in  the  wings  of  broadened  lines  can  be  an  important  absorption  mechanism.  Thus 
the  ability  to  propagate  laser  light  efficiently  through  gaseous  ablation 
products  depends  on  whether  the  products  absorb  the  laser  radiation  and  on  the 
fate  of  the  absorbed  energy. 

3.2.2  Self-Consistent  Calculation  of  Gas  Heating  by  Laser  Absorption 

In  the  previous  report  we  presented  expressions  to  calculate  the  average 
number  density  nt  in  the  gaseous  plume.  Here  we  shall  restate  these 
expressions  in  order  to  explain  how  we  have  solved  the  several  equations 
self-consistently  to  predict  the  temperature  of  the  plume.  This  solution  is  an 
approximation  where  we  consider  the  gas  in  the  region  of  one-dimensional  ( 1 -D ) 
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flow  near  the  ablating  surface.  Here  we  solve  for  the  pressure  and  temperature 
of  a  mixture  of  ideal  gases  in  the  approximation  that  the  pressure  and  tempera¬ 
ture  are  uniform  over  the  region  from  the  edge  of  the  Knudsen  layer  to  the  end 
of  the  region  of  1-D  flow,  that  the  thickness  of  the  1-D  region  is  increasing 
at  the  sonic  velocity,  and  that  the  species  concentrations  are  always  at  their 
equilibrium  values.  The  input  parameters  include  the  incident  laser  irradiance 
I0  for  the  CW  case  (I0  =  Fp/fp  for  pulsed  operation),  and  spot  diameter  D,  as 
well  as  the  elemental  composition  of  the  ablating  material,  the  absorptivity  of 
the  surface  a,  and  the  absorption  cross  sections  a  of  the  major  absorbing 
species  in  the  plume.  The  solution  is  iterated  until  the  temperature  and 
pressure  in  the  uniform  plume  are  consistent  with  the  amount  of  energy  absorbed 
in  the  plume  and  at  the  target  surface. 

The  rate  of  material  ablated  from  the  target  can  be  estimated  using  the 
heat  of  ablation  (J/g)  calculated  for  equilibrium  from  the  amount  of  energy 
absorbed  per  unit  area  as 

dm/dt  =  a  Is/H^  .  (7) 

Here  dm/dt  is  the  time  derivative  of  the  mass  removed  per  unit  area,  evaluated 
after  passing  through  the  Knudsen  layer  where  we  assumed  that  the  directed 
velocity  is  sonic  after  the  Knudsen  layer.  Thus  dm/dt  is  evaluated  from  the 
equilibrium  vapor  pressure  calculation  as  described  in  our  previous  report 
(Simons,  et  al.,  1989,  Appendix  A).  In  the  above,  a  is  the  fraction  of  the 
surface  irradiance  Is  that  is  absorbed  at  the  target.  The  density  of  material 
in  the  gas  p  (assuming  all  products  are  gaseous)  is  given  by 


58 


(8) 


p  =  (dm/dt)/ug  =  a  Is/(ug  HA)  , 

where  ug  is  the  flow  velocity,  which  we  take  to  be  sonic,  so  that 

Ug  =  (Y  Rg  T/<M>)0-5  =  9.82  x  105  (y  T/<M>)0-5  cm/s  ,  (9) 

where  Rg  is  the  gas  constant,  T  is  the  temperature  (Kelvin),  y  is  the  ratio 
of  heat  capacities  Cp/cv,  and  <M>  is  the  average  molecular  weight. 

The  fraction  fA  of  the  incident  irradiance  that  is  absorbed  by  the  a 
mixture  of  gaseous  ablation  products  is  determined  by  the  total  absorption 
cross  section  ,  the  total  number  density  nt,  and  the  path  length  L^d  as 

f A  =  Is/Io  =  exP("CTt  nt  l1D>  *  (10> 

where  I0  is  the  incident  irradiance.  The  is  the  sum  of  the  neutral  average 

absorption  cross  section  <a>  and  the  cross  sections  from  the  two  inverse 
bremsstrahlung  processes  of  electron-neutral,  with  cross  section  <Jibn,  and 
electron-ion,  with  The  average  <a>  is  determined  from  the  absorption 

cross  sections  of  the  different  species  and  the  mole  fractions  x^  as 

<cr>  =  £  cr-[  x^  ,  (11) 

and 

fft  =  [<cr>nt  +  cribn  (1  -  xe)  xe  nt2  +  aibi  xe2  nt2l 

x  [1  -  exp(-Ep/Te) ]  ,  (12) 

where  the  second  term  in  brackets  corrects  for  stimulated  emission,  Ep  is  the 
photon  energy  (1.165  eV),  and  Te  is  the  temperature  in  eV.  Here,  xe  is  the 
mole  fraction  of  electrons,  which  we  assume  to  equal  the  mole  fraction  of  ions, 
and 

nt  =  p  Na/<M>  ,  (13) 

where  NA  is  Avogadro's  number  (6.02  x  1022  molecules/mole)  and 
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(14) 


<M>  =  E  Xi  Mi  , 

where  is  the  molecular  weight  of  species  i. 

The  effective  thickness  (i.e.,  the  path  length)  of  the  one-dimensional 
region  of  the  ablation  plume  for  short  pulse  times  Tp  can  be  estimated  as 

Lid  =  ^g  *^p  *  *tp  ^  0. 55D/ Ug  ,  (15a) 

when  the  pulse  time  is  short  enough  that  the  expansion  is  one-dimensional.  For 
long  pulse  times,  the  expansion  become  three-dimensional  and  the  value  of 
depends  on  the  diameter  of  the  illuminated  target.  As  discussed  in  our 
previous  report  (Simons,  et  al.,  1989),  Bailey  found  that  the  true  ( p  L) 
product  integrated  over  the  distance  from  the  target  can  be  approximated  by  the 
p  from  Eq.  (8)  times  a  characteristic  length  given  by 

L1d  s  0.55D  ,  Tp  >  0. 55D/ug  ,  (15b) 

where  D  is  the  diameter  of  the  laser  beam  at  the  target  surface. 

In  order  to  include  the  absorption  by  inverse  bremsstrahlung  processes,  we 
need  the  electron  density,  ne,  and  values  for  the  absorption  cross  sections. 

We  use  tha  Saha  equation  to  find  the  mole  fraction  of  electrons,  xe  =  ne/nt,  as 
xe  =  [6.04  x  1021  (1  -  xe)  nt(gl/g0)  T1*5  exp(-Ipe/Te) ]°- 5/nt  ,  (16) 

For  our  calculations,  we  take  g\/gQ  =  6/9  and  Ipe  =  (11.26-1.17)  eV,  values 
that  are  appropriate  for  atomic  carbon,  although  they  do  not  change  much  when 
we  consider  the  Cn  or  C2H  species.  Note  that  we  use  an  effective  ionization 
potential  Ipe  that  is  reduced  from  Ip  by  the  value  of  one-photon  of  energy. 

This  is  to  account  for  the  rapid  ionization  expected  at  high  pressures  for  any 
species  excited  to  within  1.17  eV  of  its  ionization  potential. 

Ve  find  the  cross  section  Pibn  for  the  electron-neutral  inverse 
bremsstrahlung  process  by  interpolation  among  the  results  calculated  for 
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different  X  and  T  for  atomic  carbon  by  Geltman  (1973).  For  example,  at  5000  K, 

cribn  =  1-23  x  10~39  cm5  at  \  =  l.o  ym  and  T  =  5000  K.  We  compute  the  cross 

section  ffibi  for  the  electron-ion  inverse  bremsstrahlung  process  from  the 

Kramers  formula  times  the  Gaunt  factor  <gff>  as 

cjibi  =  ( 4/ 3 )  [2  ix/(3  m  kg  T )  J 0 * ^  [z2  e6/(b  c  m  \>3)]<gff>  , 

or 

aibi  =  1.27  x  10-37  X3  z2  <gff>/Te0-5  cm5  ,  (17) 

where  z  is  the  charge  of  the  ion,  and  the  constant  is  for  X  given  in  ym  and  the 

temperature  Te  in  eV.  For  our  case  of  initial  plasma  formation,  z  =  1.  We 

find  <gff>  by  interpolating  between  values  given  as  a  function  of  variables 
determined  by  X  and  T  at  z  =  1  (Karzas  and  Latter,  1961-2).  As  <gff>  is  on  the 

order  of  one  and  Oibi  is  around  100  o^bn ,  both  processes  contribute  when  xe  is 

greater  than  0.01  or  so. 

The  final  step  in  the  calculation  is  to  find  the  new  temperature  from  the 
increase  in  enthalpy  when  the  gas  absorbs  energy.  The  increase  in  enthalpy  of 
the  gas  from  absorption  of  a  constant  fraction  fA  of  the  incident  irradiance  I0 
over  a  time  At  is  fAI0At.  The  mass  of  gas  in  the  plume  is  a(l-fA)I0At/HA,  so 
that  the  enthalpy  change  per  unit  mass  is 
AHg  =  [fA/(l  -  fA)]HA/<*  . 

When  the  energy  in  the  directed  flow  is  subtracted,  this  becomes 

AHg  =  [fA/(l  -  fA)](HA/a)  -  ug2/2  .  (18) 

The  total  enthalpy  of  the  gas  is  then  AHg  +  HA.  We  then  use  the  results  of 
equilibrium  thermochemical  calculations  that  give  Hg  as  a  function  of  T  and  the 
pressure  p. 

When  the  procedure  finds  self-consistent  values  for  p,  T,  xj,  etc.,  then 
the  amount  of  material  ablated  is  consistent  with  the  absorption  in  the  plume, 
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and  the  temperature  in  the  plume  is  consistent  with  the  energy  absorbed  there. 
The  iteration  procedes  as  shown  in  Figure  30.  The  procedure  converges  quickly 
when  <  0.2  or  so,  especially  if  the  absorption  is  primarily  by  the  neutral 
species.  At  high  absorption  via  inverse  bremsstrahlung,  the  procedure  may  not 
converge.  Therefore,  we  somewhat  arbitrarily  declare  plasma  formation  under 
conditions  where  the  temperature  is  increasing  and  f^ib  *s  substantial. 

3.2.3  Absorption  Line  Width  Estimation 

The  importance  of  absorption  by  isolated  lines  depends  on  the  absorption 
width  of  the  lines  and  on  the  number  density  of  the  absorber.  As  we  are 
primarily  interested  in  plasma  initiation,  where  the  electron  number  density  is 
low,  we  ignore  the  line  broadening  from  the  Stark  effect.  This  leaves  three 
sources  of  line  broadening:  natural,  Doppler,  and  collisional.  The  natural 
width  is  determined  by  the  lifetime  of  the  upper  state  via  the  Heisenberg 
uncertainty  principle.  A  lifetime  of  5  x  10“^s  corresponds  to  a  width  of  0.01 
cm--*-,  which  is  small.  Thus  lifetime  broadening  is  only  important  when  the 
lifetime  is  shorter  than  10-*-®s.  The  Doppler  half  width  is  determined  by  the 
frequency  of  the  radiation  and  the  velocity,  which  is  a  function  of  the 
temperature  and  the  mass, 

yd  =  3.58  x  10-8  \>(T/M)0-5  ,  (19) 

where  \>  is  the  frequency  (in  cm-*-)  of  the  transition,  T  is  the  temperature 
(Kelvin),  and  M  is  the  gram  molecular  weight.  The  collisional  or  Lorentz 
broadening  is  a  function  of  the  collision  cross  section  and  the  number  density. 
Thus  the  Lorentz  half  width  is  given  by 
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Begin  at  Low  Irradiance 


Given:  Fp,  Tp,  I0,  D,  c^(T,p) 

Preliminary  Calculation  of  HA(<xIs),  Hg(T,p), 
and  Xi(T,p)  from  Elemental  Composition  of  Target 

Initially  Assume:  xe  =  0,  use  T  and  p  for  Target  Surface  at  I0 

[1]  Find  <M>  and  y  at  equilibrium  for  given  T,  p  from  x^(T,p)  tables. 
Evaluate 

<ct>  =  E  oj  xj  ,  (11) 

at  =  [<a>nt  +  aibn  (1  -  xe)  xe  nt2  +  cjibi  xe2  nt2] 

x  [1  -  exp(-Ep/Te) ]  .  (12) 

[2]  Find  nt  by  iteration  of 

f A  =  Is/Io  =  exP(~°t  nt  l1D>  »  <10> 

nt  =  a  IQ( 1  -  fA)NA/(<M>  ug  HA)  ,  (7)+(8)+(13) 

until  nt  consistent  with  fA. 

[3]  Update  the  inverse  bremsstrahlung  absorption  by  finding  a  new  xe  from 

xe  =  [6.04  x  1021(1  -  xe)  nt(g1/g0)  T1-5  exp(-Ipe/Te)]0-5/nt  .  (16) 

[4]  Loop  back  to  [2]  to  update  the  at  value  using  the  new  xe  and  nt  to  find 

the  new  cxt  using  cribn  and  at  the  current  T.  Stop  loop  when  change 

in  fA  is  small. 

[5]  Use  fA  and  HA  to  find  Hg. 

[6]  Find  new  T  from  Hg(T,p). 

[7]  If  AT  >  0.5  K,  loop  back  to  [1]  using  the  new  T. 

[8j  If  procedure  converged,  increment  Fp  and  go  to  [ 1 ] .  Use  most  recent 
values  as  initial  values  for  next  point. 


Figure  30. 

Iteration  Procedure  to  Self-Consistently  Solve  for  Plume 
Temperature  as  a  Result  of  Heating  by  Gas  Absorption 
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YL  =  n  p(T/300  )-m  ,  (20) 

where  p  is  the  pressure  in  atm  and  and  m  depend  on  the  collision  cross 
section  with  typical  values  ca.  0.1  cm“l  atm--*-  and  0.5  to  1.0,  respectively. 

The  number  density  above  an  ablating  target  surface  is  also  a  function  of 
temperature,  and  Figure  31  shows  combined  effects  of  the  Doppler  and 
collisional  half  widths  for  a  graphite  epoxy  target  when  m  =  0.7  and 
n  =  0.1  cm-l  atm--*-. 

Molecular  species  have  absorption  bands  that  normally  have  many  absorption 
lines  to  excite  individual  rovibrational  levels  in  the  bound,  upper  electronic 
state.  These  bands  can  also  form  a  pseudocontinuum,  and  this  is  especially 
likely  if  the  gas  temperature  is  high  so  that  many  vibrational  levels  are 
populated.  For  example,  the  C3  absorption  from  370  to  490  nm  studied  at  3200  K. 
had  lines  so  closely  spaced  that  they  could  not  be  resolved  with  a  spectrograph 
resolution  approaching  the  Doppler  width  of  a  single  line  (Brewer  and  Engelke, 
1962).  The  pseudocontinua  can  be  understood  as  the  result  of  lines  spaced  more 
closely  than  their  line  width,  and  they  can  also  occur  in  diatomic  species  when 
the  pressure  is  high  enough,  as  we  show  in  the  next  subsection. 

3.2.4  Calculation  of  C?  Phillips  Band  System  Spectrum 

The  C2  radical  has  a  number  of  electronic  transitions  that  can  be  excited 
with  photon  energies  near  the  1.165  eV  laser  photon  energy.  For  this  diatomic 
radical,  the  absorption  spectrum  in  the  IR  at  low  pressures  will  consist  of  a 
series  of  a  lines  representing  transitions  betwen  specific  states  specified  by 
their  rotational,  vibrational,  and  electronic  quantum  numbers.  Therefore,  the 
a  value  at  the  laser  wavelength  will  depend  very  strongly  on  the  laser 
wavelength,  as  the  a  will  be  large  if  the  absorption  line  happens  to  overlap 
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LINE  WIDTH 


Figure  31. 

Line  Half  Width  as  a  Function  of  Pulse  Fluence  for  a 
Graphite  Epoxy  Target  with  Tp  =  60  and  300  ns 
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the  laser  line.  As  the  pressure  increases  the  linewidth,  the  absorption  lines 
will  overlap,  and  the  dependence  of  a  on  X  will  be  reduced.  As  we  discussed 
previously  (Simons  et  al.,  1989,  Subsection  3.3.8),  the  potential  energy 
curves  of  C2  show  a  ground  state  X(*Eg+),  which  has  an  allowed  transition  with 
an  excitation  energy  of  1.025  eV  to  the  A(lnu)  state.  This  transition  is 
called  the  Phillips  band  system.  In  addition,  there  is  a  another  IR  system, 
the  Ballik-Ramsay  system,  with  transitions  from  the  lowest  triplet  states,  the 
a(3nu)  level  at  an  excitation  energy  of  0.0756  eV  relative  to  the  ground  state 
to  the  b(3Eg~)  level  at  an  excitation  energy  (relative  to  the  a(^nu)  state)  of 
0.698  eV.  There  is  another  triplet  system  between  the  c(^Eu+)  level  with  a 
total  excitation  energy  of  1.657  eV,  which  should  have  allowed  transitions  to 
the  d(3jlg)  level  located  0.821  eV  higher.  The  possibility  of  several  elec¬ 
tronic  transitions  from  states  at  energies  of  0.076  to  1.66  eV  gives  a  reason¬ 
able  expectation  that  C2  absorption  can  be  important,  especially  when  the 
pressure  is  high. 

In  this  subsection  we  describe  calculations  of  the  C2  absorption  spectrum 
for  one  of  the  IR  transitions,  the  Phillips  band  system  (X  ^Eg+  to  A  ^nu).  The 
calculation  used  a  modification  of  the  first  part  of  the  PSI  spectral-fitting 
code.  The  input  parameters  for  the  code  are  the  spectral  constants  for  the 
upper  and  lower  states,  the  Einstein  A  coefficients  for  the  transitions  between 
the  various  vibrational  levels  in  the  two  electronic  states,  the  temperature, 
and  a  specification  of  the  line  shape  to  be  used.  The  calculations  are  for  the 
12 C2  isotopic  molecule.  We  calculated  the  Einstein  A  coefficients  for  transi¬ 
tions  between  the  vibrational  quantum  levels  of  v'=0  to  15  in  the  upper  elec¬ 
tronic  state  going  to  vibrational  levels  of  v"=0  to  19  in  the  lower  state. 
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This  calculation  requires  values  for  the  Franck-Condon  factors  qV'V"  between 
the  different  vibrational  wave  functions,  as  well  as  the  electronic  transition 
moment  Re,  because 

A  =  T<Re>2[64ltV(3h\3)]  [  (2- 6q f  A,  +  A»  )/(2- Sq,  A' )  ]  ,  (21) 

where  h  is  the  Planck  constant  and  the  second  term  in  brackets  is  2  for  Z-Z 
transitions,  and  1  for  other  cases,  as  A'  and  A"  are  the  angular  momenta  of  two 
electronic  states  (Whiting,  et  al.,  1980;  Schadee,  1978). 

The  calculation  of  the  r-centroids  and  Franck-Condon  factors  necessary  to 
evaluate  the  A  coefficients  used  a  modification  of  the  PSI  code  that  Marinelli 
used  to  calculate  A  factors  for  N2,  which  has  been  described  previously 
(Marinelli  and  Piper,  1985).  Rydberg-Klein-Rees  (RKR)  potentials  were  first 
calculated  for  both  electronic  states  using  the  spectroscopic  constants  of 
Huber  and  Herzberg  (1979),  and  the  approach  of  Tellinghuisen  (1974).  This  code 
was  modified  to  fit  the  portion  of  the  potential  between  the  internuclear 
separations  r  at  the  classical  turning  points  of  the  lowest  vibrational  state 
with  a  modified  Morse  potential.  The  potential  used  was  that  of  Hulbert  and 
Hirschfelder ,  as  described  in  Herzberg  (1950),  which  has  the  form 

U(r)  =  De[(l  -  e-fr)2  +c(33xV2$c(l  +  bf3x)]  ,  (22) 

with  x  =  r  -  re 

|3  =  [2n2  cy/(Deh)]0*5  afe 

c  =  1  -  [1  +  0^/(6  Be2)]/ (f5re) 

b  =  2  -  {(7/12)  -  [(5/4)  +  5agC^/(12Bg2) 

+  5og2(^2/(144Be4)  -  2afexe/(3Be)]/(f32re2)}/c  .  (23) 

Note  the  corrected  sign  for  the  second  term  of  the  expression  for  b  in 
Eq.  (23),  which  is  discussed  in  Hulbert  and  Hirschfelder  (1961).  In  the  above, 
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re  is  the  equilibrium  internuclear  separation,  y  is  the  reduced  mass,  and  Be, 
«g,  and  <%xe  are  spectrosopic  constants.  We  did  the  fit  of  the  negative 
and  positive  x  values  separately  by  adjusting  the  dissociation  energy  De  so 
that  U  at  the  x-value  of  the  classical  turning  point  agrees  with  the  potential 
from  the  RKR  calculation. 

The  numerical  eigenfunctions  resulting  from  this  procedure  were  then 
evaluated  using  the  Numerov-Cooley  procedure  (Eccles  and  Malik,  1981)  to  solve 
the  radial  schrodinger  equation.  The  overlap  integrals  in  the  r-centroid  and 
Franck-Condon  factor  calculations  were  evaluated  using  simpson's  rule.  The 
results,  which  cover  the  range  of  v'  and  v"  from  0  to  19,  are  more  extensive 
than  other  calculations  (e.g.,  Cooper  and  Nicholls,  1976).  We  take  the  E<Re>^ 
value  to  be  0.38(eao)^,  where  e  is  the  electronic  charge  and  ac  is  the  Bohr 
radius,  or  Re  =  1.567  Debye.  Cooper  and  Nicholls  (1976)  determined  that 
E<Re/(eaQ)>2  =  0.36  and  0.40  for  the  Av  =  0  and  +2  sequences,  respectively, 
while  Brewer  and  Hagen  (1979)  quote  two  other  sources  as  0.48  and  0.50  for  the 
Av  =  0  sequence.  Our  calculated  qV'V"  values  agree  rather  well  with  those  of 
Cooper  and  Nicholls  in  the  region  of  overlap,  e.g.  for  the  64  transitions  they 
reported,  40  agree  within  1  percent,  only  11  have  differences  >  4  percent,  and 
only  6  of  these  are  >  10  percent. 

Using  the  calculated  A  coefficients  in  the  PSI  spectral  synthesis  code, 
along  with  the  spectral  constants  and  the  temperature,  gives  line  strengths  for 
the  various  rotational  lines  of  each  vibrational  transition.  Each  contribution 
was  summed  over  the  line  shape,  which  is  assumed  to  be  a  Voigt  profile.  This 
line  shape  is  the  combination  of  the  effects  of  Doppler  and  pressure  broaden¬ 
ing.  The  code  calculates  the  power  emitted  in  Watts  per  steradian  per  micron 
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value  of  the  rotational  temperature  Tr.  The  power  from  each  vibrational  level 
is  summed  after  weighting  by  the  fractional  population  at  Tr  (i.e.,  the  elec¬ 
tronic  and  vibrational  temperatures  are  assumed  to  be  the  same  as  the  rota¬ 
tional  temperature).  The  total  power  spectrum  is  then  smeared  over  the  finite 
linewidth  of  the  Voigt  profile.  The  result  is  the  spectral  radiance  Wx,  the 
power  radiated  per  wavelength  interval  at  each  wavelength  by  each  molecule  of 
C2  at  an  equilibrium  temperature  of  Tr,  The  absorption  cross  section  per 
molecule  is  then  obtained  by  noting  that  the  absorptivity  is  related  to  the 
emissivity,  and  the  emissivity  for  one  molecule  of  C2  is  the  ratio  of  the 
calculated  spectral  radiance  to  that  radiated  by  a  blackbody  at  the  same 
temperature  and  wavelength,  which  is 

Wb5  =  3.745  x  104X~5{exp[1.4388  x  104/(XT)  ]-l}~1  ,  (24) 

when  X  is  expressed  in  ym,  T  is  in  Kelvin,  and  Wbb  is  in  W/(cm2  ym).  That  is, 
the  absorption  is  related  to  the  emissivity  t  by 

e  =  Wx/Vb5  •  <25> 

and 

e  =  1  -  exp  (-a  n  L)  =  a  ,  (26) 

where  the  second  equality  in  Eq.  (26)  considers  that  exp(-cr  n  L)  is  given  by 
1  —  cr  when  we  consider  one  molecule  of  C2. 

Figure  32  shows  the  spectrum  from  1.01  to  1.10  ym  at  T  =  5000  K,  which  is 
a  forest  of  closely  spaced  lines  when  the  line  width  is  small.  This  figure  is 
for  a  pressure  of  1  atm,  and  the  line  shape  was  calculated  using  Eqs.  (19)  and 
(20)  with  m  =  0.7  and  =  0.19  cm~l  atm-l.  This  value  of  y\  results  from 

assuming  that  the  appropriate  collision  cross  section  for  the  line  broadening 
is  the  same  as  the  cross  section  from  viscosity  studies.  We  take  the  cross 
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Figure  32. 


Absorption  Spectrum  Calculated  for  C2  X(-*-Eg+)  to  A(-*-IIu)  from  1.01  to 
1.10  pm.  The  calculation  is  for  T  =  5000  K  and  p  =  1  atm,  and 

amav  =  5.06  X  10-16  cm2 


section  as  (ji/4)(ctq  +  )  >  where  cr^  and  are  the  radii  given  by  the 

Lennard-Jones  parameters  for  C  and  C2  (Svehla,  1962).  The  spectrum  includes 
contributions  from  v"=  0  to  19,  going  to  v'=  1  to  15  in  the  upper  state. 

Figure  33  shows  the  region  from  1.06  to  1.07  pm  on  an  expanded  scale.  This 
figure  shows  how  increasing  the  pressure  at  T  =  6000  K  changes  the  spectrum 
from  one  of  isolated  lines  at  1  atm,  to  lines  that  are  slightly  overlapping  at 
10  atm,  to  lines  that  are  strongly  overlapped  at  100  atm,  with  the  structure 
nearly  gone  at  300  atm.  Figure  34  shows  how  increasing  the  pressure  to  100  atm 
changes  Figure  32.  Thus,  the  modulation  of  the  spectrum  depends  on  the 
pressure.  The  structure  of  the  spectrum  does  not  depend  very  strongly  on  the 
temperature,  as  shown  in  Figure  35.  This  shows  the  spectra  at  p  =  10  atm  for 
T  =  4000  and  7000  K,  which  have  nearly  the  same  structure  in  this  region.  The 
magnitude  of  a  at  the  particular  value  of  X  =  1.064  ym  is  weakly  dependent  on 
temperature,  and  it  does  depend  on  the  line  width.  This  dependence  is  strong 
when  the  line  half  width  is  <  0.7  cm~l,  as  shown  in  Figure  36.  Here  the 
different  symbols  represent  T  values  from  4000  to  7000  K.  The  line  width  is 
primarily  determined  by  the  pressure,  and  30  atm  corresponds  the  half  widths  of 
0.7  to  1  cm~l.  Figure  37  shows  the  data  of  Figure  36  plotted  versus  pressure. 
In  summary,  this  calculation  for  a  (k-  1.064  ym)  gives  values  that  are  large 
and  approximately  independent  of  pressure  for  p  >  30  atm,  and  a  increases  with 
pressure  at  lower  pressures. 

3.2.5  Estimation  of  C3  Absorption  Spectrum 

The  C3  radical  is  a  major  species  in  the  equilibrium  vapor  above  carbon- 
rich  targets,  as  shown  later  in  Figure  42.  In  our  previous  report  (Simons,  et 
al.,  1989)  we  showed  an  extrapolation  of  the  data  of  Brewer  and  Engelke  (1962) 
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Figure  33. 

Effect  of  Pressure  on  the  Absorption  Spectrum  Calculated  for  C2  X(-*-Eg+)  to 
A(illu)  from  1.06  to  1.07  ym  at  T  =  6000  K.  In  parts  (a)  through  (dj,  the 
pressure  is  1,  10,  100,  and  300  atm,  +  yl  is  0.0836,  0.294,  2.39, 
and  7.06  cm-*,  and  amax  is  33.2,  13.1,  3.09,  and 
2.10  x  10“17  cm^,  respectively. 
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Figure  35. 

Effect  of  Temperature  on  the  Absorption  Spectrum  Calculated  for  C2  X(1E„+)  to 
A(1nu)  from  1.06  to  1.07  pm  at  p  =  10  atm.  In  parts  (a)  and  (b),  the 
temperature  is  4000  and  7000  K,  Yp  +  yt  is  0.359  and  0.275  cm~  ,  and 
is  1.19  and  1.28  x  10~16  cm2,  respectively. 
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Figure  36. 

Effect  of  Line  Half  Width  (yp  +  y^)  on  the  Absorption  Cross  Section  Calculated 
at  1.064  ym  for  C2  X(^Eg+)  to  A(^nu).  The  diamonds,  squares,  ovals,  and 
triangles  repesent  T  =  4000,  5000,  6000,  and  7000  K,  respectively. 
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Figure  37. 


Effect  of  Pressure  on  the  Absorption  Cross  Section  Calculated  at  1.064  pm  for 
C2  X(1Eg+)  to  A(!nu).  The  diamonds,  squares,  ovals,  and  triangles  repesent 
T  =  4000,  5000,  6000,  and  7000  K,  respectively. 
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for  absorption  to  the  lowest  excited  singlet  state  A(^IIU),  whose  excitation 
energy  from  the  ground  state,  v"=0  level,  is  24,676  cm-l  (3.06  eV,  or 
0.405  ym) .  This  comparison  failed  to  include  the  absorption  data  of  Cooper  and 
Jones  (1979),  who  show  that  the  transition  extends  out  to  at  least  0.53  ym  when 
the  temperature  is  3238  K.  Figure  38  compares  the  data  of  Cooper  and  Jones  at 
three  temperatures  with  the  Brewer  and  Engelke  data  that  we  converted  to  a 
values  by  using  the  vapor  pressure  data  of  Leider  and  coworkers  (1973).  The 
comparison  shows  good  agreement  between  the  two  experiments  at  the  maximum  of 
the  absorption,  and  the  Cooper  and  Jones  data  falls  off  faster  at  low  energies. 
Note  that  the  Cooper  and  Jones  data  at  different  temperatures  show  a  higher 
maximum  at  lower  temperatures,  and  higher  low-energy  wings  at  higher  T.  This 
is  consistent  with  the  expectation  that  the  absorption  will  spread  over  a 
greater  energy  range  as  T  increases  and  that  the  total  oscillator  strength  for 
the  transition  is  conserved.  Thus  a  constant  area  under  the  o-E  curve  is  spread 
over  a  larger  energy  range  as  T  increases. 

Figure  38  also  shows  extrapolations  of  the  data  to  1.064  ym  (1.165  eV) 
under  the  assumptions  that  the  transition  moment  Re  and  the  Franck-Condon 
factors  are  all  constant,  so  that  the  dependence  on  photon  energy  E  is  simply 
<*E)  =  c(E0)  (E/E0)  exp[-(E0  -  E)/(kB  T)]  ,  (27) 

where  Eq  is  the  photon  energy  at  the  last  measured  value,  kg  is  the  Boltzmann 
constant,  and  T  is  the  temperature.  The  figure  shows  symbols  plotted  at  the 
lowest  reported  energies,  and  the  same  symbol  at  1.165  eV  represents  the 
extrapolation  based  on  the  lowest-energy  experimental  point  at  that 
temperature.  The  figure  shows  that  the  measured  fall-offs  have  similar  slopes 
from  the  maxima  to  the  lower  energy  points  as  the  slope  between  the  last 
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measured  point  and  1.165  eV,  which  is  given  by  Eq.  (27).  This  means  that  the 
Franck-Condon  factors  are  roughly  constant  over  the  energy  range  from  about  2.4 
to  3.1  eV.  However,  we  have  no  evidence  that  the  factors  remain  constant  for 
the  remainder  of  the  extrapolation  to  1.165  eV.  As  discussed  previously 
(Simons  et  al. ,  1989),  it  is  quite  possible  that  the  vertical  transitions 
between  the  potential  curves  have  a  larger  separation  than  1.165  eV,  which 
would  mean  that  the  absorption  would  have  to  be  a  non-vertical  transition  with 
a  small  Franck-Condon  factor.  Thus  the  extrapolations  should  be  considered  to 
give  upper  limits  to  the  a  value  at  1.064  pm. 

Figure  38  also  shows  an  estimate  for  the  absorption  resulting  from 
transitions  between  the  two  lowest  triplet  states.  The  first  electronically 
excited  state  of  C3  is  a  triplet  state,  a(3nu),  which  is  observed  at  2.12  eV  in 
a  Ne  matrix  (Jacox,  1988),  a  value  that  would  correspond  to  about  2.10  eV  in 
the  gas.  This  state  can  be  thermally  populated  at  high  temperatures  in  the 
plume.  For  example,  the  degeneracy  times  the  Boltzmann  factor  is  0.013,  0.044, 
and  0.10  at  T  =  4000,  5000,  and  6000  K,  respectively.  The  next  excited  triplet 
state  is  the  3jig  level  that  has  an  excitation  energy  relative  to  the  lowest 
triplet  predicted  to  be  0.71  to  0.81  eV  (Romelt,  et  al.,  1978),  0.93  eV 

(Williams,  1975),  or  0.74  eV  (Jones,  1985).  Thus,  there  is  an  allowed 

transition  predicted  to  be  accessible  from  the  lowest  triplet  state.  The 
oscillator  strengths  for  the  triplet  transitions  are  not  available,  but  we 

would  expect  the  maximum  a  values  to  be  in  the  range  of  10“^-^  to  10“^  cm“3. 

In  the  figure,  we  have  estimated  a  at  T  =  4000,  5000,  and  6000  K  by  assuming 
that  the  square  of  the  transition  moment  Re^  and  the  Franck-Condon  factors  are 
the  same  as  they  are  for  the  X(^Ig+)  to  A(*nu)  transition.  What  we  plot  is  the 
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Figure  38. 

Absorption  Spectrum  of  C3.  The  high  energy  region  shows  data  at  3200  K 
(Brewer  and  Engelke,  1962),  adjusted  to  use  vapor  pressures  of  Leider  et  al 
(1973),  and  data  at  T  =  3238,  3720,  and  4012  K  (Cooper  and  Jones,  1979). 
These  measurements  are  extrapolated  to  give  the  points  shown  at  1.165  eV 
(see  text).  Also  shown  is  an  estimate  of  the  low-energy  triplet  transition 
a(3nu)  to  b(-^IIg)  for  T  =  4000,  5000,  and  6000  K.  The  top  curve  shows 
the  corresponding  a  values  if  the  a(^nu)  state  were  the  ground  state. 
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Cooper  and  Jones  (1979)  oqj  data  at  3720  K  adjusted  using  Eq.  (27)  and  the 
thermal  population  of  the  a(^nu)  state  to  give 

cr3(E)  =  aCJ(E  +  0.74)  [E/(E  +  0.74)]g3  exp(-2 . 12/Te)  ,  (28) 

where  g3  is  the  degeneracy  factor  of  six.  The  upper  curve  shows  what  a  would 
be  if  the  a(^nu)  state  were  the  ground  state  (i.e.,  the  g3  exp(-2. 12/Te)  term 
is  set  to  unity)).  This  shows  a  maximum  o  of  lO-^-^  cm^,  which  is  a  reasonable 
value.  It  is  possible  that  the  the  a  to  b  transition  will  have  sharp  struc¬ 
ture,  so  that  a  at  the  laser  line  could  be  a  small  value  between  the  transi¬ 
tions.  The  fact  that  the  lowest  singlet  transition  is  a  pseudocontinuum  at 
T  =  3200  K,  argues  that  the  triplet  excitation  will  also  be  pseudocontinua  at 
high  T,  although  the  lowest  frequency  vibration  for  the  a  state  is  predicted  to 
be  larger  than  that  for  the  ground  state,  which  will  spread  the  line  spacing. 

At  higher  pressures,  the  absorption  will  become  a  pseudocontinuum  in  the  same 
way  as  we  saw  above  for  C2-  Figure  39  shows  data  from  Figure  38  in  a  plot 
versus  wavelength.  Here  the  X  to  A  <r(1.064  pm)  values  using  Eq.  (27)  to 
extrapolate  the  data  from  the  two  experimental  groups  are  shown  at  T  =  4000, 
5000,  and  6000  K,  with  the  experimental  point  used  in  the  extrapolation  shown 
with  the  same  symbol  as  the  extrapolated  points. 

3.2.6  Absorption  Spectra  of  Other  Plume  Species 

In  graphite  epoxy  and  Kevlar  epoxy  targets,  the  C2H  radical  is  an 
important  species  in  absorbing  1.06-pm  radiation.  It  has  a  low-energy 
electronic  transition  from  the  X(^E+)  to  the  A(^I1)  state,  with  a  transition 
energy  predicted  to  be  around  0.5  eV  (Shih,  et  al.,  1977,  1979),  and  the  band 
at  0.446  eV  has  been  tentatively  identified  as  the  origin  (Yan,  et  al. ,  1987). 
This  is  a  triatomic  radical  with  a  low-frequency  bend,  which  leads  to  the 
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Figure  39. 

Absorption  Spectrum  of  C3.  The  high  energy  region  shows  data  from  Brewer 
and  Engelke  (1962)  at  3200  K  and  the  T  =  3238  and  4012  K.  data  of 
Cooper  and  Jones  (1979).  The  ovals  and  the  triangles  represent 
extrapolations  from  the  measured  data  (at  the  X  shown  by  the  point) 
to  the  points  shown  at  1.064  pm  for  T  =  4000,  5000,  and  6000  K 
(see  text).  Also  shown  is  an  estimate  of  the  lowest  triplet 
transition,  a(^nu)  to  b(3llg)  for  T  =  4000,  5000,  and  6000  K. 
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expectation  that  absorption  will  be  a  pseudocontinuum  at  higher  temperatures, 
epecially  when  the  pressure  is  large  enough  to  substantially  broaden  the  lines. 
Reimers  and  coworkers  (1985)  have  calculated  the  absorption  spectrum  at  3000  K. 
The  maximum  absorption  of  <7  =  3  x  10~19  cm2  near  1.064  pm. 

The  C3H  and  C4H  radicals  are  also  important  species  in  the  plume  from 
graphite  epoxy  targets.  It  is  reasonable  that  these  will  also  have  sizable  a 
values  at  1.06  urn,  as  these  will  have  at  least  one  low  frequency  vibration, 
which  implies  the  absorption  spectrum  will  contain  pseudocontinua  at  high 
temperatures.  In  addition,  there  are  several  diatomic  species  that  are 
important  in  ablation  plumes,  and  the  absorption  spectra  for  these  should  be 
calculated  in  the  manner  described  above  for  C2 . 

3.2.7  Equilibrium  Calculation  of  Significant  Species  in  the  Ablation  Plume 

At  the  high  gas  pressures  that  exist  near  the  surface  of  a  target  that  is 
rapidly  ablating,  it  is  reasonable  to  assume  that  the  reaction  rates  are  fast 
enough  that  the  species  are  at  chemical  equilibrium,  epecially  when  we  consider 
the  high  temperatures  in  the  region.  The  equilibrium  values  we  use  are  calcu¬ 
lated  in  the  same  way  as  earlier  calculations  at  PSI  (Simons,  et  al.,  1989, 
Appendix  A),  except  that  we  have  included  more  species.  In  particular,  we  have 
included  ionic  species,  so  that  the  enthalpy  continues  to  increase  with  T,  even 
when  there  are  essentially  no  molecular  species  left.  This  is  shown  in 
Figure  40,  which  gives  Hg  as  a  function  of  T  for  pressures  from  10”^  to  10^  atm 
for  a  graphite  epoxy  target.  The  relative  species  concentrations  for  this 
material  are  shown  as  a  function  of  the  temperature  in  Figures  41  and  42  for 
pressures  of  10  and  1000  atm,  respectively.  These  show  that  molecular  species 
are  important  up  to  higher  temperatures  when  the  pressure  is  higher.  At  high 
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Figure  40. 

Enthalpy  per  Gram  of  the  Equilibrium  Gas  Mixture  That  Has  the 
Elemental  Composition  of  Graphite  Epoxy  at  Different  Pressures 

as  Functions  of  Temperature 
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Figure  41. 


Mole  Fractions  of  Species  in  the  Equilibrium  Gas  Mixture  That  Has 
Elemental  Composition  of  Graphite  Epoxy  at  p  =  10  atm  as 

Functions  of  Temperature 


the 
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MOLE  FRACTION 


Figure  42. 

Mole  Fractions  of  Species  in  the  Equilibrium  Gas  Mixture  That  Has  the 
Elemental  Composition  of  Graphite  Epoxy  at  p  =  1000  atm  as 

Functions  of  Temperature 
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pressures  the  important  molecular  species  are  C4H,  C3,  C2H,  and  C2 ;  also,  CO, 
C3H,  C4,  C5,  Cg,  C7,  CH,  H2,  C2H+,  and  C2H2  are  present  with  x-^  >  0.01  for  some 
temperature  region.  The  equilibrium  calculations  depend  on  the  thermochemical 
values  used  as  input.  The  results  presented  here  use  values  that  are  primarily 
from  the  JANAF  tables  (Chase  et  al.,  1985),  as  discussed  previously  (Simons, 
et  al.,  1989,  Appendix  A).  The  major  features  of  the  results  are  similar  when 
we  use  a  thermochemical  data  set  developed  for  combustion  studies  at  Sandia 
National  Laboratory  (Miller  and  Bowman,  1988).  This  data  set  is  a  revision  of 
one  developed  by  NASA  workers  (see  Burcat,  1984). 
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